High pt resonances as a possibility to explore hot and dense nuclear matter 
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One of the fundamental objectives of experiments with ultrarelativistic heavy ions is to explore 
strongly interacting matter at high density and high temperature. In this investigation we study in 
particular the information which can be obtained by analyzing baryonic and mesonic resonances. 
The decay products of these resonances carry information on the resonances properties at the space 
time point of their decay. We especially investigate the percentage of reconstructable resonances as 
a function of density for heavy ion collisions in the energy range between Eiab = 30 AGeV and ^/s 
— 200 AGeV, the energy domain between the future FAIR facility and the present RHIC collider. 

PACS numbers: 24.10.Lx,25.75.-q,25.75.Dw 



The experimental analysis of heavy ion reactions us- 
ing resonances has been applied for several years from 
low energy [ll, |3| through intermediate [1, 0| to high en- 
ergy heavy ion collisions 0, S 0] • general one distin- 
guishes between leptonic and hadronic decay channels. 
While the hadronic decay channels have the advantage 
of larger branching ratios, the leptonic decays have the 
advantage that the decay particles do not undergo final 
state interactions. Thus it is worthwhile to work out the 
differences and the advantages and disadvantages of the 
two approaches which will be discussed in the following. 

The present Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven and the upcoming Facility for Antipro- 
ton and Ion Research (FAIR, for a recent status on the 
project we refer to @) provide an excellent research envi- 
ronment for probing resonances in matter. At the RHIC 
experiments it has been observed Q that less resonances 
are measured than expected from statistical model calcu- 
lations ■ Stable hadrons however follow the prediction 
of this model. This suggests the conclusion that after 
chemical freeze-out, when the chemical composition of 
the final state is determined, hadrons still undergo col- 
lisions and therefore some of the resonances cannot be 
identified by the invariant mass of the decay products. 

At FAIR the leptonic as well as the hadronic decay 
channel can be explored. While the leptonic channel is 
usuall y re garded as the 'cleaner' channel recent calcula- 
tions have shown that the dilcpton channel might 
not probe the dense phase as it was expected before. 
In light of this new development it is worthwhile to eval- 
uate the density-profile and the space-time-evolution of 
resonances which can be reconstructed in the hadronic 
decay channels. Although those channels suffer from the 
drawback of final state interaction of the decay prod- 
ucts, their large branching ratios might make them better 
suited for the investigation of the high density phase of 
heavy ion collisions compared to leptonic decay channels. 

For our calculations we utilize the UrQMD(v2.3) 



model, a non-equilibrium transport approach, which re- 
lics on the covariant Boltzmann equation. All cross sec- 
tions arc calculated by the principle of detailed balance 
and the additive quark model or arc fitted to available 
data. UrQMD does not include any explicit in-medium 
modifications for vector mesons or effects to describe 
the restoration of chiral symmetry. The model allows 
to study the full space time evolution of all hadrons, 
resonances and their decay products in hadron-hadron 
or nucleus-nucleus collisions. This permits to explore 
the emission patterns of resonances in detail and to 
gain insight into their origins and decay channels. For 
previous studies of resonances within this model see 
tni, m, El, Q E [H, [I3- For further details about 
the UrQMD model the reader is referred to pl.[T9j. 

Experimentally, the reconstruction of resonances is 
challenging. One often applied technique is to recon- 
struct the invariant mass spectrum for single events. 
Then, an invariant mass distribution of mixed events is 
generated (here, the particle pairs are uncorrelated by 
definition). The mixed event distribution is substracted 
from the invariant mass spectrum of the single (corre- 
lated) events. As a result one obtains the mass distri- 
butions and yields (after all experimental corrections) 
of the resonances by fitting the resulting distribution 
with a suitable function (usually a Brcit-Wigncr function 
peaked around the pole mass of the respective resonance). 

If the resonance spectral function changes in the 
hadronic medium this is in principle visible in the differ- 
ence spectrum between true and mixed events. However, 
if a daughter particle (re-)scattcrs before reaching the 
detector the signal for the experimental reconstruction is 
blurred or even lost. Especially for strongly interacting 
decay products this effect can be sizeable. It is therefore 
difficult to judge whether a deviation from an expected 
Breit-Wigner distribution is due to an initial deformation 
or an increase of the initial width or due to the momen- 
tum dependence of the rescattcring cross section of the 
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daughter particles. 

What makes this analysis even tougher is the fact that 
the resonances decay over a wide range of densities and 
therefore only an average value is measured. If this aver- 
age value is dominated by resonance decays at low density 
the information from the high density phase is blurred 
and may offer only a limited view on the high density 
phase of the heavy ion collision. 

UrQMD offers a different technique for the extraction 
of resonances which we apply here. We follow the indi- 
vidual decay products of each decaying resonance (the 
daughter particles). If the daughter particles do not 
rescatter in the further evolution of the system, the res- 
onance is counted as "reconstructablc" . The advantage 
of this method is that it allows to trace back the origin 
of each individual resonance to study their spatial and 
temporal emission pattern. Because UrQMD follows the 
space time evolution of all particles it is possible to link 
production and decay point of each individual resonance. 
This method also allows to explore the reconstruction ef- 
ficiency in different decay branches. 

In order to calculate at which density the resonance 
decays we have to determine the baryonic density. The 
baryon density is calculated locally at the position of the 
resonance in the rest frame of the baryon current (Eckart 
frame) as pB = j'^ with ~ (ps, 0). Details on the cal- 
culation of the baryon density are discussed in In 
all figures we present the density in units of ground state 
density, where a value of 0.16 1/fm'^ is assumed. In the 
following we discuss the density dependence of the prob- 
ability that a resonance can be reconstructed. Naively, 
one would expect that the higher the densities the more 
the rescattering effect becomes dominant. Therefore it is 
unlikely that a resonance which decays at high density is 
reconstructablc. The view on the low density zone is ex- 
pected to remain unblurred but is less interesting because 
it resembles that observed in elementary collisions. 

Depicted in Fig. [1] top (bottom) is the probability 
that a baryon resonance - shown arc A, S* (1385) and 
A* (1520) baryon resonances (p,a;, and $ mesons) - 
which was produced at a certain density can be recon- 
structed experimentally. One observes a clear peak at 
very low density and a steady decrease towards higher 
density. This means that resonances that are produced at 
rather low density have a high probability to be detected 
and as the density increases the chance to reconstruct the 
resonances decreases. This is nothing unexpected. How- 
ever, this trend stops at roughly 2 po- At higher densities 
the chance to reconstruct a resonance saturates or even 
increases slightly again. This increase, which we discuss 
later in detail, is caused by resonances which picked up 
very high transverse momenta and leave the interaction 
zone quickly. This results in a decay in a region with 
less hadronic activity and a higher chance to be recon- 
structed. 

Whereas the form of the curves is qualitatively similar 
for the different hadrons the absolute value of the fraction 
of reconstructablc resonances is rather different. It can 
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FIG. 1: (Color online) Fraction of reconstructable baryon res- 
onances (top) and meson resonances (bottom) as a function 
of baryon density at the point of production. 



be understood in terms of lifetimes of the resonances and 
in terms of the rescattering cross sections of the decay 
products. 

Due to the large cross section of pions in nuclear matter 
(usually undergoing iV + 7r^Aor7r + 7r^p reactions) 
the probability to detect a high density A resonance or 
a p meson is rather small compared to the probability to 
detect a high density $ meson, since the $ meson itself 
has a small cross section in nuclear matter and a long 
lifetime of ~ 40 fm/c and the hadronic decay products 
(mostly kaons and antikaons) have a smaller cross sec- 
tions when compared to the pions from the decay of a 
p meson. Similarly, the long lifetime of the A increases 
their possibility to be reconstructed. As mentioned ear- 
lier, the saturation or slight increase of the reconstruction 
probability as a function of density has its origin in the 
possibility that resonances with a large pt can escape 
quickly from the reaction zone which is rather small ini- 
tially. 

Fig. [2] shows for various resonances the probabil- 
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FIG. 2: (Color online) Probability distribution of baryon den- 
sity at the production vertex for various reconstructable reso- 
nances in central (b < 3.4 fm) Au-|-Au collisions at 30 AGeV 
(top figure) and 200 AGeV (bottom figure) as a function of 
baryon density. One observes that most resonances which can 
be reconstructed in the hadronic decay channel originate from 
low baryon density. 



ity that an experimentally reconstructable resonance has 
been created at a density p. The integral over all densities 
is normalized to unity. One observes that most of those 
resonances are produced at very low densities, which is 
especially true for the mesonic resonances. 

Reconstructable baryon resonances stem from slightly 
higher baryon densities, however most are still produced 
at rather low densities (with a peak at roughly 0.1 
ground state density). So the detection of resonances 
produced at densities above ground state densities using 
hadronic decay channels seems not too encouraging. 
However, as we discuss next, a loophole might exist. 

Let us illustrate this further with an example which is 
representative for all investigated particles. 

Fig. [3] depicts the average transverse momentum of 



A resonances as a function of baryon density. Lines 
show reconstructable resonances, symbols show all de- 
cayed resonances. The striking feature is the different av- 
erage transverse momentum between all resonances and 
those which are reconstructable. The higher the average 
transverse momentum, the larger is the chance that the 
resonance can be reconstructed. The < pT > of recon- 
structable A resonances is about 200 MeV higher than 
for all A resonances. Resonances with a large pr can 
leave the high density zone rather fast and move with a 
velocity of about < px > /"m- outwards. 

Another interesting feature in Fig. [3] is the differ- 
ence between the 7^=200 AGeV and Eiab=iO AGeV 
curves. While the Eiab~50 AGeV data shows a de- 
crease of < pt > as a function of the baryon density, the 
^=200 AGeV data show an increase. At ^=200 AGeV 
the initial collisions (which happen at high baryon den- 
sity) are more energetic and give the particles a high 
transverse momentum, subsequent rescattering decreases 
Pt- For the iSjafe =30 AGeV collisions the situation is op- 
posite. Initially the particle pr is small and the rescat- 
tering increases the px due to transverse expansion. 
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FIG. 3: (Color online) Average transverse momentum of re- 
constructable (line) or all (symbol) A baryons as a function 
of baryon density for two different energies. 

Fig. m shows the pr dependence of the reconstruction 
probability in detail. It shows the transverse momentum 
spectra for all (full symbols) and reconstructable reso- 
nances (open symbols) . The numbers stated in the three 
shaded areas {px < 1 GeV, 1 GeV < pr < 2 GcV, pr > 
2 GeV ) are the percentages of reconstructable resonances 
created at a density higher than 2po. One observes that 
at low transverse momentum the percentage of recon- 
structable resonances is low and increases when going to 
higher transverse momenta, i.e. that with increasing px 
the chance to reconstruct a resonance produced at high 
baryon density increases. 

In conclusion, we have discussed that the view on the 
high density zone may not be as restricted as usually 
assumed when analyzing hadronic resonances. 
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We argued that resonances detected with high trans- 
verse momentum are sensitive to higher densities. It will 
be interesting to explore if the properties of these reso- 
nances are different from the bulk emitted at low densi- 
ties. The exploration of high pT resonances might there- 
fore open a new keyhole at the upcoming CBM exper- 
iment at FAIR or the critRHIC program to gain infor- 
mation on the high density zone and to observe eventual 
changes of resonance properties in the medium. 



FIG. 4: (Color online) Transverse momentum spectra for all 
and reconstructable resonances for central (b<3.4 fm) Au-|-Au 
collision at 30 AGeV beam energy. Full circles depict the 
spectrum for all decayed resonances (included in the analysis 
are A, A, E baryons, as well as p,u), K"^ and u mesons), open 
circles for reconstructable resonances. The numbers indicate 
the percentage of reconstructable resonances stemming from 
density region with p/po > 2. 
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